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ABSTRACT 



The star Cyg OB2 No. 9 is a well-known non-thermal radio emitter. Recent theoretical work suggests that all such O-stars should be in a 
binary or a multiple system. However, there is no spectroscopic evidence of a binary component. Re-analysis of radio observations from the 
VLA of this system over 25 years has revealed that the non-thermal emission varies with a period of 2.35 ± 0.02 yr. This is interpreted as a 
strong suggestion of a binary system, with the non-thermal emission arising in a wind-collision region. We derived some preliminary orbital 
parameters for this putative binary and revised the mass-loss rate of the primary star downward from previous estimates. 
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1. Introduction 

The Cygnus OB2 association has been the focus of many 
massive-star studies. Observations show that it contains more 
than 2600 OB stars, of which ~ 120 are O stars (Knodlseder 
2000). It also harbours some of the hottest and most luminous 
stars known in our Galaxy. Several of these stars exhibit signifi- 
cant variations in their X-ray and/or radio continuum emission. 

One of these stars is Cyg OB2 No. 9 (VI Cyg 9, Schulte 9, 
MT 431), which is classified as an 05If star (Massey & 
Thompson ll991l >. It was one of the first O stars to be identified 
as a non-thermal radio source (Abbott et al. 1984). It shows the 
classical properties of non-thermal radio emission, such as a 
high brightness temperature (Phillips & Titus 1 1990] ) and, also, 
a radio spectral inde?(Q with a characteristic negative value. 

The non-thermal emission comes from relativistic elec- 
trons, believed to be accelerated in shocks. Although shocks 
are ubiquitous in O-star winds due to the instability of the ra- 
diative driving mechanism (e.g. Owocki & Rybicki 1984), the- 
oretical studies show that this embedded-shock model cannot 
be the source of the non-thermal emission, and that the accel- 
erating shocks arise where the two stellar winds collide (Van 



Send offprint requests to: S. Van Loo 
e-mail: svenvl@ast . leeds .ac.uk 

1 The radio spectral index a describes the power-law behaviour of 
the flux: F v oc V oc A~ a . For thermal emission a ss +0.6. 



Loo et al. 120061 ). This conclusion necessarily implies that all 
non-thermal emitting O stars are in a binary system. A similar 
correlation between non-thermal radio emission and binarity is 
already well established for WR stars (Dougherty & Williams 
I2000t which are the evolutionary descendants of O stars. 



At present, there is no spectroscopic evidence that Cyg 
OB2 No. 9 is a binary. De Becker ( 120051 ) and Kiminki et 
al. (|2006) measured the radial velocity on a number of occa- 
sions. Neither study found significant variation on time-scales 
of a few months to a year. In the absence of spectroscopic evi- 
dence of binarity, we need to look to other indicators. Possible 
techniques include high-spatial resolution radio observations 
that reveal the characteristic shape of a wind-collision region 
as observed in WR 140 (Dougherty et al. 120531 ) and WR 147 
(Williams et al. 11 9901 ), or searching for radio variability consis- 
tent with the orbital motion of a binary (e.g. WR140: Williams 
etal. H990V 



In this paper we demonstrate that the radio emission in Cyg 
OB2 No. 9 has a clearly periodic behaviour. Radio data ob- 
tained with the Very Large Array (VLA) since 1980 were re- 
reduced, revealing a 2. 35-yr period. It is argued that this is evi- 
dence of a binary system, with the non-thermal radio emission 
arising in a wind-collision region. 
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Fig. 1. The observed 6-cm fluxes of Cyg OB2 No. 9 as a func- 
tion of time. 



2. Data and radio light curve 

We collected archive data for Cyg OB 2 No. 9 from the NRAO 
Very Large Arraycl archive. Many of these observations have 
not been published previously. To avoid introducing system- 
atic effects in the data reduction, we re-reduced all observations 
consistently. Details of the reduction are given in Appendix A. 
The measured fluxes and their error bars (or upper limits) are 
listed in Table A. 1 . 

Using the observations listed in Appendix lAl we can look 
for variability at radio wavelengths. In Fig. Q] we see that the 
6 cm fluxes of Cyg OB2 No. 9 are highly variable with maxima 
found at e.g. 24 May 1980 and 3 July 1982, and minima at e.g. 9 
May 1983 and 19 September 1985. These extrema give a rough 
estimate of a period, that is ;/any period is present, of the order 
of 2 yr. However, due to the problem of aliasing, we cannot 
exclude any divisor of 2 yr. 

To find periodicity in the data, we apply the string-length 
technique described by Dworetsky j 1983b to the different 
wavelength observations. The period is found by minimising 
the distance between the flux observations in an orbital phase 
diagram, i.e. the string length. The string-length method is par- 
ticularly useful for a small number of observations with a large 
separation in time. Furthermore it does not require a specific 
shape for the light curve. However, this method does not work 
when multiple periods are present in the data (e.g. Bourguignon 
et al. 120071 1. Also, this technique neglects the effect of errors 
on the observation. We, therefore, modified the method by re- 
placing the flux difference , i.e. {f\ - f^) 1 , by the mean flux- 
difference, or (/i - fi) 2 + {cr 2 + cr 2 ,) (where cr is the error on the 
flux) when calculating the distance between the observations. 

We find a minimum in the string length between 2.3 and 
2.4 yr in both the 3.6 and 6 cm fluxes (see Fig. [2j. However, 
there is no periodicity found in the 2 and 20 cm fluxes. This is 
not surprising as the 2 cm fluxes have a high degree of noise 
making a periodicity search difficult. At 20 cm, we do not find 




2 The National Radio Astronomy Observatory is a facility of the 
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Fig. 2. The string length as a function of period for different 
wavelength observations of Cyg OB2 No. 9: 6 cm (solid), 3.6 
cm (dashed), 2 cm (dotted). We do not plot the 20 cm fluxes 
as it is similar to the 2 cm line. The minimum representing the 
most-likely period is found near 2.35 yr. 



a period because most observations have a similar flux level 
or a high upper limit (see Fig. [3}. Using the method described 
by Fernie (I1989I ), we can estimate the true period and its un- 
certainty. We find P = 2.34 ± 0.03 yr for the 6 cm fluxes and 
P = 2.37 + 0.03 yr for the 3.6 cm fluxes. A weighted mean and 
variance then gives P = 2.355 + 0.015 yr. This fits well with our 
crude first estimate and is confirmed by other period-finding 
methods such as a discrete Fourier transform (Scargle[l982|l. 

To assess the goodness of the period, Fig. [3] shows the 2, 
3.6, 6 and 20 cm radio fluxes of Cyg OB 2 No. 9 folded with 
the 2. 35-yr period. Phase 0.0 has been arbitrarily set at the low- 
est 6 cm upper limit (JD 2443762.18). The consistency in the 
light curve shows that these reduction techniques do not intro- 
duce systematic effects, though they may be responsible for a 
few minor effects. There are some low fluxes around phase 0.68 
at 6 cm and a somewhat higher flux at phase 0.27 that we at- 
tribute to uncertainty in the absolute flux calibration (see notes 
in Table[AD- 

The large number of 6 cm observations around phase 0.23- 
0.28 (see inset to Fig. 01 could in principle be useful to search 
for variations on shorter timescales. There does not seem to 
be significant variability however. A linear regression for this 
subset shows that, except for one observation which has a de- 
viation of 3cr, all observations can be fitted. Furthermore, the 
flux increase is consistent with the gradient seen in the main 
figure. 

3. Discussion 

3.1. Binarity 

The 2.35 yr period derived for the non-thermal radio emission 
of Cyg OB 2 No. 9 is strongly suggestive of a binary system 
with non-thermal radio emission arising in a wind-collision re- 
gion. The intrinsic synchrotron radiation of the wind-collision 
region changes as a function of the binary separation, while 



12 
10 

8 

6 

4 

2 


-( 
12 

10 
8 
6 
4 

2 


-( 
12 

10 

8 
6 
4 

2 


-( 
12 

10 

8 
6 
4 

2 


-( 



Van Loo et al.: Non-thermal radio emission from O-type stars 



T 



I 



I 



o.o 



0.2 



0.4 



0.6 



0.8 



1.0 



- I 



eI 



JeI 

I 



o.o 



0.2 



0.4 




0.6 



13 



I 



0.8 



1.0 



I 



0.4 



0.6 



0.8 



1.0 



1. 



I 



J 

I 



I 



.2 



0.0 



0.2 0.4 0.6 

phase in 2.35-yr cycle 



0.8 



1.0 



, 3.6, 6 and 20 cm radio fluxes of Cyg OB 2 No. 9, folded with a 2.35-yr period. Phase 0.0 has been arbitrarily se 
:m upper limit (JD 2443762. 18). 



4 



Van Loo et al.: Non-thermal radio emission from O-type stars 



free-free absorption removes some of the intrinsic emission, 
depending on the position in the wind (i.e. optical depth along 
the line-of-sight to the observer). Other mechanisms such as 
inverse-Compton cooling and the Razin effect are also depen- 
dent on the binary separation (e.g. Pittard et al. 120061 1. These 
mechanisms therefore produce a radio light curve that is modu- 
lated by the binary orbit. For WR 140, the best studied example 
of a colliding-wind binary, the radio light curve shows good re- 
peatability from one orbital period to another (White & Becker 
[19951 

The binarity of Cyg OB2 No. 9 could in principle show up 
in other observational diagnostics. We review a number of them 
below. 

Optical observations The semi-major axis, a, of the binary 
can be estimated from 

(a[AU]f = M[M ](P[yr]) 2 , (1) 

where M is the mass of the binary system. A lower limit on the 
mass of Cyg OB2 No. 9 can be estimated on the basis of the 
Eddington limit 

KeL * <1. 
4nGM,c 

Using log(LJL bo ,) = 6.24 (Herrero et al. [19991 we find M, > 
4OM0. For hot, massive stars, however, the Eddington parame- 
ter is often within a factor of two below unity, so that the stellar 
mass is more likely to be M, > 80Mq. This agrees well with 
observational estimates of the evolutionary mass, i.e. ~ IOOMq 
(Herrero et al. 1 19991 >. Although we cannot say anything about 
the actual mass of the binary, it is unlikely to exceed « 140Mq, 
the mass of the most massive hot-star binary known, WR20a 
(Rauw et al. 120041 ). For a binary mass between 80 and 140 
Mq, we find a semi-major axis of 7.6 AU 55 a 9.2 AU from 
Eq. CQ. 

The separation between the binary components lies within 
the (1 +e)a range. This means that, even for large eccentricities, 
the maximum separation is % 18 AU. With a distance to the Cyg 
OB2 association of 1.82 kpc (Bica et al. 120031) . the components 
have a maximum separation on the plane of only 0.01 arcsec. 
A Hubble FGS (Project: FGS 10612; PI: D. Gies[| observation 
of Cyg OB No. 9 in April 2006 looked for a binary companion 
with luminosity difference less than 4 magnitudes within a 0.01 
- 1 arcsec region around the primary. Given the estimated sep- 
aration, it is unlikely that a companion can be detected using 
FGS. 

Radial velocity measurements Two independent studies of 
radial velocity variations by De Becker (12005b and Kiminki et 
al. ((2006) did not find any significant changes in radial veloc- 
ity. However, by combining the measurements of these stud- 
ies, we find that a constant radial velocity can be rejected on 
a 97 % confidence level. Thus, Cyg OB2 No. 9 is a possible 
binary by the criterion of Kiminki et al. ( 120061 ). Furthermore, 

3 URL: http://www.stsci.edu/observing/phase2-public/10612.pro 
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eccentricity 

Fig. 4. Contour lines for sin ij Vl - e 2 , where the solid line 
corresponds to a 70 + 70 Mq binary and the dashed to a 
100 + 40 Mq binary. The permitted range of eccentricity and 
inclination lies within the contour lines. 

re-examination of spectroscopic archive data shows a possi- 
ble double-lined spectrum (Kiminki, priv. comm.). Additional 
spectroscopic observations are needed to confirm the prelimi- 
nary spectroscopic evidence of binarity. 

With the derived constraints on the binary separation and 
the period, the expected amplitude of the radial velocity varia- 
tions can be calculated. The observed variation is likely to be 
due to the motion of the primary. For a Keplerian orbit, the 
radial velocity amplitude is given by 

2na p sin/ 2na M s sin; 

" ~ p vT^? = P M P + M s VT^ 2 ' 

where a p and a are the semi-major axis of the primary star and 
the binary system respectively, M p and M s the mass of the pri- 
mary and secondary star respectively, i the inclination angle 
and e the eccentricity. For a 100 + 40 Mq binary, a p = 0.28a 
which gives K p = 32 sin ij Vl - e 2 kms 4 . For a 70 + 70 Mq 
binary, one finds K p = 55 sin// Vl - e 1 km s _1 . This can then 
be compared to the observed radial velocity amplitude. By 
combining both data sets of De Becker ( 120051 ) and Kiminki 
et al. (120061 . we find a marginally significant variation with a 
radial velocity amplitude K p between 27 and 69 km s _1 . 

We find that this comparison allows a first estimate of the 
orbital parameters of the system, through the constraint on the 
function sin// Vl - e 2 . Figure [4] shows the permitted values of 
the eccentricity and the inclination for different binary mass 
ratios. Due to the broad range of possible radial velocity am- 
plitudes, we only find a weak constraint on the inclination and 
eccentricity. However, the derived value for the radial velocity 
variation is subject to considerable uncertainty: the sampling 
of the orbit is sparse, which makes it likely that the maximum 
variation was missed. More recent observations suggest that 
the radial velocity amplitude tends to the upper value of our 
adopted range (Naze, pers. comm.). This favours a high eccen- 
tricity orbit (e > 0.6) viewed at a substantial inclination, e.g. 
/ > 60°. 
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X-ray emission A non-thermal radio emitting binary should 
produce both thermal and non-thermal X-rays. The non- 
thermal component is produced by inverse-Compton cooling 
of relativistic electrons in the strong stellar ultraviolet radia- 
tion field. The thermal component is produced by the hot gas 
(~ 10 6 K) produced in the wind collision. For a long period bi- 
nary such as Cyg OB2 No. 9, this thermal component is hard, 
and dominates the non-thermal contribution in the 0.5 - 10 
keV range (De Becker [20051 . Depending on orbital parame- 
ters, some modulation of the X-ray emission with the orbital 
period can be expected. However, the X-ray emission of Cyg 
OB2 No. 9 remains relatively stable for ~ 15 years (Waldron 
et al. 11998b . However, it should be noted that the X-ray ob- 
servations are sparsely distributed throughout the orbital phase 
space, and variations due to orbital changes might have been 
missed. 

By fitting the X-ray data with a two-temperature model, 
Rauw et al. (12005b find that the higher temperature of the two, 
i.e. 12xl0 6 K, agrees well with the temperatures expected in 
a wind-collision binary (Stevens et al. 119921 . Interestingly, the 
observed X-ray emission from No. 9 is smaller than the canon- 
ical Lx/L bo i = 10~ 7 (Waldron et al. |1998b . Although some scat- 
ter around this value can be expected, binaries are often much 
brighter in X-rays than single stars (Pollock ll9871 l. In this re- 
gard, Cyg OB2 No. 9 does not have a typical X-ray luminosity 
for a binary. 

Emission line variability Another possible indicator of collid- 
ing winds is the shape and variability of the emission lines. In 
the colliding-wind binary Cyg OB2 No. 5, the He II ^4686 line 
shows orbital phase variations, but with a phase offset com- 
pared to the photospheric lines. This shows that at least part 
of this line is formed in the colliding-wind region (Rauw et 
al. |1999l ). Similar effects could be present in Cyg OB2 No. 9. 

High-spatial resolution radio observations Recent high- 
resolution radio observations of Cyg OB2 No. 9 with the 
VLBA reveal a bow-shaped emission region (Dougherty et al., 
in prep). This is consistent with emission arising in a wind- 
collision region, and is very reminiscent of the emission dis- 
tribution observed in proto-typical colliding-wind systems like 
WR 140 (Dougherty et al. 120051 1. This is the strongest direct ev- 
idence to date that the non-thermal radio emission in Cyg OB2 
No. 9 arises due to a colliding-wind binary. 

3.2. Spectral index 

Fig. [5] shows the fluxes at chosen orbital phases, as a function 
of the available wavelengths. Globally, the radio spectrum be- 
comes shallower as the orbit progresses from radio minimum 
(phase ip ~ 0) to radio maximum (<p « 0.65 Near radio min- 
imum (<p = 0.97), the spectral index between 2 and 6 cm is 
q-26 > 0.26, suggesting an optically-thick spectrum with sub- 
stantial attenuation, most likely due to free-free absorption. 

4 Recall that a positive slope on the figure corresponds to a negative 
spectral index, and vice versa. 
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Fig. 5. The radio spectrum at chosen orbital phases. A thermal 
radio spectrum and a pure synchrotron spectrum are plotted at 
arbitrary flux levels, for the sake of reference. The numbers 
refer to the orbital phase of the observation. 



One may surmise that the minimum occurs near periastron, 
where the wind-collision region is buried in the stellar wind, 
though we readily admit this is not the only possible config- 
uration for a minimum. The shape of the spectrum near radio 
minimum suggests thermal emission, but does not exclude a 
modest non-thermal contribution. Between phase 0.1 and 0.3 
the spectral index a-26 switches from positive to negative. The 
spectral index between 6 and 20 cm is always negative. At radio 
maximum the spectrum has a turnover from negative to posi- 
tive spectral index as one goes from short to long wavelengths. 
This is consistent with the attenuation expected from free-free 
absorption, which is greater at larger wavelengths. The spec- 
tral index at maximum is q-26 = -0.24 ± 0.09, which differs 
from the -0.5 index expected for synchrotron emission arising 
from strong-shock acceleration. This suggests that, even at ra- 
dio maximum, the observed flux is a combination of optically 
thick and optically thin emission. 

3.3. Mass-loss rate 

We can derive an upper limit for the mass-loss rate from the 
radio minimum, where we can expect the contribution of non- 
thermal emission to be smallest. Adopting the stellar parame- 
ters given in Van Loo et al. (|2004) and taking the 0.3 mJy upper 
limit at 6 cm (1 1 June 1990), we find M < 0.9 x 10~ 5 M© yr" 1 . 
The 0.8 mJy detection at 2 cm gives M < 1.1 X 10~ 5 M© yr -1 . 
These values are upper limits, as the observed flux at radio 
minimum may still have a non-thermal contribution. They are 
lower than the M(radio) = 1.9 xlO" 5 M© yr~' given by Abbott 
et al. (11984b . partly because of the lower observed flux, and 
partly because of improved stellar and wind parameters. 

The mass-loss rate originally derived from Ha 
is 1.77xl0~ 5 Moyr 1 (Leitherer [19881 . which yields 
1.47xl0~ 5 M© yr~' with the parameters given in Van Loo 
et al. (12004b . A recent study by Puis et al. (12006b shows 
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that, for O star supergiants with Ha in emission, M(radio) 
~ 0.4 - 0.5 M(Ha). This is believed to be caused by a radial 
decrease in the amount of clumping, as is also predicted 
theoretically (Runacres & Owocki 120051 1. We find for Cyg 
OB2 No. 9 that M(radio) < 0.61M(Ha), which is roughly 
consistent with the results of Puis et al. ( 120061 1. 

3.4. Energy budget of the non-thermal radio emission 

The radio synchrotron emission is maintained by the kinetic 
power of the stellar winds, but only part of that is dissipated 
through the wind-collision shocks, and only a fraction of that 
power is dissipated in relativistic electrons. In addition, radio 
synchrotron emission accounts for a small fraction of the total 
non-thermal energy, with the bulk of that energy radiated at 
high energies. 

There is considerable uncertainty how much of the wind- 
collision power is converted into radio synchrotron luminos- 
ity. Following Chen & White ( 11994t and Pittard & Dougherty 
(2006 ) we suggest the following order-of-magnitude estimate 
of 

L syn ss 10 4 LNTe ~ 5 x 10 7 Lwc, 

where L^j e is the power in non-thermal electrons and Lwc is 
the kinetic power of the winds that enter the wind-collision re- 
gion. 

The kinetic power Lwc depends strongly on the momentum 
ratio of the winds. For equal winds, the kinetic power in the 
wind-collision region is roughly half of the total kinetic power 
of the winds, whereas for a system with a much smaller mo- 
mentum ratio e.g WR 140, the kinetic power in the winds is 
only 1 % of the total wind power due to the small solid angle 
presented by the wind-collision region to the dominant wind 
(Pittard & Dougherty [2006b . For Cyg OB2 No. 9, using the 
stellar wind parameters, we estimate that L syn « 10 31 erg s _1 
if the wind momenta are assumed to be equal, and L sync « 
10 29 erg s- 1 if the companion has a much weaker wind. The 
observed radio synchrotron luminosity ~ 10 30 erg s _1 , which 
is likely a lower limit on the intrinsic synchrotron emission due 
to attenuation of emission by absorption and/or cooling. In sup- 
port of these numbers, the model of Eichler & Usov ([1993 ) 
gives a similar value. Given all the uncertainties in the above 
discussion, these numbers are all closely compatible, support- 
ing a wind-collision origin for the non-thermal emission. 

4. Conclusions 

The most important result of this paper is that a period of 2.35 
yr was found in the radio emission from Cyg OB2 No. 9. This is 
highly suggestive of a binary system in which the non-thermal 
emission, modulated by the orbit, arises in a wind-collision re- 
gion. The binary status of Cyg OB2 No. 9 can currently not 
be confirmed by spectroscopic, optical or X-ray observations. 
However, some of the observations support the presence of a 
binary e.g. the high temperature derived from the X-ray spec- 
trum and the marginally significant radial velocity variations. 
Additional observations are required to substantiate the claim 
of a binary using these techniques. However, there is already 



compelling evidence of a wind-collision region based on high- 
resolution radio observations (Dougherty et al., in prep). 

From the radio minimum we derive an upper limit for the 
mass-loss rate of 0.9xl0~ 5 Mq yr~ ! . This is substantially lower 
than previous estimates. As the radio formation region is likely 
to be less affected by clumping than the formation region of 
Ha, we suggest our upper limit is the most reliable estimate to 
date. 
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Appendix A: Data reduction 

The data reduction was accomplished using the Astronomical 
Image Processing System (AIPS), developed by the NRAO, 
following the standard procedure of antenna gain calibration, 
absolute flux calibration, imaging, and deconvolution. These 
have been described in detail in Paper I. 

For some of the observing programmes that we examined, 
no flux calibrator was observed and we used calibration obser- 
vations from observing programmes at closely the same time. 
For programmes VP5 1 and AL372 at 20cm no such observa- 
tions were available and we resorted to the value listed in the 
VLA calibrator manual making the fluxes from these obser- 
vations less reliable. As a consistency check on the derived 
flux scale, we compared fluxes for the most commonly used 
phase-reference source, J2007+404, with those obtained at 2, 
3.6, and 6 cm at the University of Michigan Radio Astronomy 
Observatory (UMRAO). Excellent agreement was found for all 
observations except AC1 16 (1985-01-13) at 2cm. For this spe- 
cific case, the flux scale used was derived from interpolation of 
the MRAO data. 

Antenna gain calibration was established in the standard 
manner of self-calibration of the phase-reference source and 
transferring these solutions to the target field-of-view. These 
solutions were then improved further by one round of phase- 
only calibration of the target field-of-view. In some cases, 
smoothly varying phase solutions were not found using the ad- 
ditional target self-calibration, and so we adopted the phase- 
referenced solutions. These cases where target self-calibration 
was applied are noted in Table I A. fl 

Images of the field-of-view were derived by filtering out 
the lowest spatial frequency data. This eliminates large scale 
"background" emission that is particularly prevalent at 20cm 
in the more compact 'C and 'D' configuration observations. 
This should have no effect on the flux estimates of Cyg OB2 
No. 9 since this source is unresolved by the VLA in the highest 
resolution 'A configuration at all wavelengths examined. 

After image deconvolution using AIPS task CLEAN, the 
flux of Cyg OB2 No. 9 was established by fitting elliptical 
Gaussians to the source image. The resulting fluxes are listed 
in Table lA. ll Correction for primary beam attenuation has been 
made for observations where Cyg OB2 No. 9 is not at the field 
centre. Some care has been taken to assess the uncertainty in 
the derived fluxes. The uncertainties given in Tabl dA.ll include 
not only the root-mean-square (RMS) noise in the map, but 
also an estimate of the uncertainty in the absolute flux cali- 
bration process (5% at 0.7 and 2cm, and 2% at 3.6, 6, 20 and 
90 cm) and an estimate of systematic errors. The systematic er- 
rors were evaluated using a jackknife technique, described in 
Paper I. In those observations where Cyg OB2 No. 9 was not 
detected, an upper limit is quoted that is 3 times the uncertainty 
as derived above. 

A comparison between the fluxes derived here and those 
in the literature show good agreement overall, with sometimes 
slightly higher values reported here, perhaps due to the use 
of phase self-calibration on the target field. There are a few 



notable exceptions to this agreement. In the AC42 data at 2 
and 6cm, there is weak (~ 3cr) source at the location of Cyg 
OB2 No. 9, but other stronger "sources" nearby do not repeat 
in other 2 and 6cm observations. Hence we assign a 3<x upper 
limit (contrary to Bieging et al. |19891 >.In the BIGN data, we find 
a 4cr detection, where Bieging et al. provide a 3cr upper limit. 
We do not detect Cyg OB2 No. 9 in VP91 at 20cm as reported 
by Phillips & Titus (1 19901 1. This is a D-configuration observa- 
tion with considerable large scale emission around the position 
of Cyg OB2 No. 9. It is not possible to totally eliminate this 
emission by removing the low spatial frequency data since it 
leaves insufficient data for successful imaging, but examination 
of 20cm data from higher-resolution configurations shows that 
this is background emission unrelated to Cyg OB2 No. 9. We 
speculate that Phillips & Titus may have mis-identified some of 
this emission as Cyg OB2 No. 9. In the AR328 0.7cm data, we 
find a weak "source" at the position of No. 9, but again there are 
brighter nearby sources that do not appear at other frequencies. 
Hence we have adopted an upper limit contrary the detection 
reported by Contreras et al. (1996). On the other hand, we do 
detect No. 9 at both 3.6 and 6cm from AR328. 

Lastly, we note a number of reporting errors. Bieging et 
al. (11989b list incorrectly some observations dated 1982 Feb 9. 
Also, Waldron et al ( 119981 1 give their 1991 observation as being 
at 6cm, rather than 3.6cm, and their 1992 data as taken on Jun 
24, rather than Jan 24. 



5 URL: http://www.vla.nrao.edu/astro/calib/manual 
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Table A.l. Reduction of the VLA data. Column (1) gives the programme name, (2) the date of the observation, (3) the source on which th 
observation was centred, (4) the phase calibrator name (J2000 coordinates), (5) the phase calibrator flux (Jy) and (6) distance to the observatio 
centre (degrees), (7) the integration time (in minutes) on the source, (8) the number of antennas that gave a usable signal, (9) the configuratio 
the VLA was in at the time of the observation, (10) the measured flux (in mjy) and (11) refers to the notes. Many of the VLA observations wer 
made in two sidebands, each of which has a bandwidth of 50 MHz; the exceptions are noted in column (11). Upper limits are 3 x the RMS 
Numbers between brackets in the notes column give references for those observations that have already been published in the literature. 
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(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


progr. 


date 


ctr. 


phase calibrator 




int. 


no. 


config. 


flux 


notes 








name 


flux 


dist. 


time 


ants. 




(mjy) 




0.7 cm 






















AR277 


1 t~\r\ a n a 1 n 

1994-04- 17 


9 


2007+404 


1.18 ±0.03 


4.9 


32 


9 


A 


< 6. 


(4) 


AR328 


1995-04-27 


9 


2007+404 


1.27 ± 0.01 


4.9 


23 


10 


D 


< 4. 


(4) 


1.3 cm 






















AW515 
2 cm 


1 f\f\f\ r\f no 

1999-06-08 


9 


2007+404 


1.46 + 0.02 


4.9 


40 


19 


AD 


3.5 ± 0.6 




BIGN 


1 no i 1 f\ *t f 

1981-10-16 


9 


2007+404 


4.8 ± 0.6 


4.9 


29 


27 


C 


4.1 ± 1.2 


1x50,(2) 


BECK 


1 noo r\c o 

1982-05-22 


9 


2007+404 


6.1 ± 1.3 


4.9 


13 


27 


A 


< 12. 


1x50 


BECK 


1 noo n*~7 ai 

1982-07-03 


9 


2007+404 


4.50 ± 0.08 


4.9 


9 


25 


A 


5.5 ± 0.9 


1x50,(1) 


BECK 


1982-12-21 


9 


2023+318 


3.18 ±0.04 


9.6 


14 


26 


D 


6.2 ± 0.8 


1x50 


AC42 


1 n o o r\c n n 

1983-05-09 


9 


2007+404 


3.76 ± 0.06 


4.9 


21 


27 


C 


< 1.3 


1x50,(2) 


AB228 


1 r\m no o c 

1983-08-25 


9 


2007+404 


6.9 ± 0.7 


4.9 


19 


26 


A 


3.1 ± 0.9 




AA29 


1 no A f\ A f\ A 

1984-04-04 


9 


2007+404 


3.53 ± 0.08 


4.9 


22 


27 


C 


4.8 ± 0.3 


S,(2) 


AO 16 


\ f\Ct A 11 O ^ 

l 984- ll -27 


9 


2007+404 


3.75 ± 0.06 


4.9 


38 


25 


A 


6.5 ± 0.4 


(2) 
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1984-12-21 


9 


2007+404 


3.96 ± 0.07 


4.9 


10 


27 


A 


6.5 ± 0.5 


S,(2) 
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1 no c A1 11 

1985-01-13 


9 
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3.3 ±0.2 


4.9 


28 


25 


A 
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4.9 
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20 
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3.10 ±0.04 


4.9 


14 


22 


A 


4.3 ± 0.4 




AH395 
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4.9 


29 


22 


A 


5.2 ± 0.4 
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9 
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5 
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7.2 ± 0.2 
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9 


2007+404 
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4.9 


1 2 


24 
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9 
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1.87 ±0.01 
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99 


23 
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9 


2007+404 


3.20 ± 0.06 


4.9 


19 


24 


A 


2.6 ±0.1 


S 
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9 


2007+404 


3.33 ± 0.07 


4.9 


23 


27 


A 


2.5 ±0.1 


S 
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1993-01-29 


9 


2007+404 
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4.9 


24 


27 
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2.8 ±0.1 


s 
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9 


2007+404 
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4.9 


6 


27 
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1993-02-14 


9 
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3.41 ± 0.03 


4.9 


7 


27 
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3.3 ± 0.1 


s 
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9 


2007+404 


3.8 ±0.1 
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14 


24 


B 
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3.28 ± 0.01 


4.9 


25 


27 
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O 

y 
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A 


1 1 
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n 

LJ 
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9 
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5 


27 


D 


6.7 ± 0.3 


s 


AW515 
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9 
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4.9 


10 


19 
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S.l 
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9 
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2.20 ± 0.01 


4.9 
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27 


B 


1.32 ±0.05 
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2000-06-26 


9 


2007+404 


2.40 ± 0.01 


4.9 


20 


26 


DnC 


3.2 ±0.2 




BD110 


2005-11-15 


9 


2007+404 


2.29 ± 0.02 


4.9 


21 


22 


D 


6.3 ± 0.2 


s 


AB1195 


2006-01-26 


8 


2007+404 


2.37 ± 0.05 


4.9 


31 


16 


AD 


<8. 




AB1210 
6 cm 


2006-05-24 


8 


2007+404 


2.21 ± 0.02 


4.9 


27 


21 


BnA 


7.1 ±0.8 




CHUR 


1980-05-23 


9 


2007+404 


4.73 ± 0.03 


4.9 


49 


22 


AD 


7.6 ±0.2 


1x50,(2) 


CHUR 


1980-05-24 


9 


2007+404 


4.58 ± 0.03 


4.9 


50 


23 


AD 


8.0 ± 0.7 


1x50,(2) 


BIGN 


1981-10-16 


9 


2007+404 


4.16 ±0.04 


4.9 


17 


27 


C 


6.2 ± 0.4 


lx50,S,(2) 


BECK 


1982-03-27 


9 


2007+404 


4.74 ± 0.01 


4.9 


211 


27 


A 


7.8 ± 0.2 


1x50,(1) 


BECK 


1982-05-22 


9 


2007+404 


4.59 ± 0.03 


4.9 


5 


24 


A 


7.9 ± 0.4 


lx50,S 


BECK 


1982-07-03 


9 


2007+404 


4.95 ± 0.02 


4.9 


5 


26 


A 


8.1 ±0.4 


lx50,S,(l) 


BIEG 


1982-08-26 


9 


2007+404 


4.4 ±0.1 


4.9 


9 


24 


B 


6.5 ± 0.4 


1x50,(2) 
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Table A.l. continued. 



ni 


C71 
y-i 


("31 
y j i 


(4) 


(5) 


l£,\ 

\y>) 


\ i ) 


yoi 


C01 
y 7) 


C101 


CI 1 1 


progr. 






phase calibrator 




int 




config 


flux 










name 


flux 


dist. 




ants 








BECK 


1087-17-71 




7 


2023+318 


2.86 ± 0.01 


ft 


o 


24 


j) 


S 1 4- 

j . 1 in u.7 


1 x^n s 

1 AJU,0 


AP47 


1 083 ns no 

i 70 j uj u7 


q 

7 


2007+404 


4.46 ± 0.02 


4 

4-. y 


1 
1 u 


if, 
zu 




1 7 

^ 1 .z 


1 vsn (i\ 

1 AJU,^Z ) 


AP47 


1 083 n8 9? 
i yo J uo -zz 


q 

7 


2007+404 


4.85 ± 0.01 


d 

4-. y 


1 4 


9S 

Z J 


A 


17 4-01 

1 . / IC u. 1 


(1 1 


AB228 


1 083-n8-78 

170J UO iO 


q 

y 


2007+404 


4.90 ± 0.03 


d 


1 8 
1 o 


27 




1 8 4- n 1 

1.0 in u. 1 




AR7S7 


1 083 i n 3n 

1 70 J 1U JU 


a 

7 


2007+404 


4.84 + 0.02 


4 

-4. 7 


q 

7 




A 


9 c _i_ n 9 
Z.O It u.z 


e 

.7 


Anio 


1 084. n3 D4. 


q 

7 


2007+404 


4.27 ± 0.01 


d 
-4. 7 


1 3 
i j 


°7 

z. / 


PnR 


< 7 4- n 3 
j . / m u. J 




A A78 


1 084. rn no 

1 70t -U J U7 


Q 
o 


2007+404 


4.33 ± 0.01 


d 

-4. s 


73 

Z J 


zu 


PnR 


^ 4- n 4 
j .7 m xj.H 




A A70 


1 084 04 04 

1 70t - U i T-U i + 


q 

7 


2007+404 


4.21 ± 0.02 


d 


1 i 
1 1 


77 
z / 


c 


f, s 4- n 9 

U.J It u.z 


J.lzJ 


v 1VU7 


1 084 1 f) 1 7 

1 70t -1U- 1 / 


q 

7 


2202+422 


2.72 ± 0.01 


1U. / 


1 1 
i J 


zu 


1 7 


7 8 4- n 8 
/ .0 It u.o 


c 


API 1 ft 

11U 


1 084 1 1 97 

1 70t 1 1 Z I 


q 

7 


2007+404 


4.50 ± 0.01 


4 

4-. 7 


1 

1 7 


ZJ 


A 


8 S 4- n 9 
O.J It u.z 




API 1 ft 

11U 


1 084 17 91 


q 

7 


2007+404 


4.63 ± 0.01 


d 
4-. 7 


q 

7 


°7 


A 


7 7 4- n 9 
/ . / It u.z 


Jilzj 


AF1 n? 


1 08s ns no 

i yo j uj U7 


+ 


2022+616 


2.43 ± 0.01 


9n 4 

ZU.4 


q 

7 


.zU 


D 


s 4- n 9 

J .7 It U.Z 


c 


A A47 


1 08^ nft 1 9 
i yo j uu i z 


q 

7 


2007+404 


3.87 ± 0.01 


4. 

4-. J 
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D 
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q 
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3.53 ± 0.01 


d 
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C71 
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1 os^ no 1 o 
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q 

7 
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3.56 ± 0.01 


d 

4-. J 


1 l 
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Z J 
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VM1 1 S 

V 1V1 1 1 J 
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q 

7 
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1 
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n 

z 


A 
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V. U.J 
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q 
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S 7 
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ZJ 
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1 ' 
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AB671 
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q 
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4 
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1 

1 7 


24 
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.7 


ARft7 1 


1 003 m 94 

1 77 JU 1 


q 

7 


2007+404 


3.08 ± 0.03 


d 
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77 
z / 


A 
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c 
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Table A.l. continued. 
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Notes: 



Column (3): source on which the observation was centred is mostly indicated by Schulte dl956t Cyg OB 2 numbers. Others indicators are: 
'98' is in between No. 8A and No. 9; '+' is listed in the observing log as 2031+41 1, 'X3' as Cyg X-3, 'XA as Cyg X-3A, '23' as C23, 
'80' as F80XF, '21' as 2032+41, '91' as 20315+41291, '32' as 20330+41132, '20' as 2033+40, 'J2' as 2033+4118.2 and 'SE' as TEVSE. 



1x50 


(or a similar code): bandwidth different from standard 2x50 MHz 


C 


flux calibration based on MRAO data (see text) 


1 


stripes are visible on this image, complicating the flux measurement 


PH 


flux calibration on phase calibrator 


S 


selfcalibration applied 


X 


flux calibration is very uncertain; check on flux of phase calibrator suggests factor of 2 


(1) 


White & Becker (TW^l 


(2) 


Bieging et al. JT9891 


(3) 


Phillips & Titus d 199011 


(4) 


Contreras et al. (11996) 


(5) 


Waldronetal. (11998k 



List of Objects 

'Cyg OB2 No. 9' on page[0 
'WR 140' on page[l] 
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'WR 147' on pagejl] 

'WR20a' on page@ 

'Cyg OB2 No. 5' on paged 



